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ABSTRACT 

We describe simulations of the response of a gaseous disc to an active spiral potential. 
The potential is derived from an N-body calculation and leads to a multi-armed time- 
evolving pattern. The gas forms long spiral arms typical of grand design galaxies, 
although the spiral pattern is asymmetric. The primary difference from a grand-design 
spiral galaxy, which has a consistent 2/4-armed pattern, is that instead of passing 
through the spiral arms, gas generally falls into a developing potential minimum and 
is released only when the local minimum dissolves. In this case, the densest gas is 
coincident with the spiral potential, rather than offset as in the grand-design spirals. 
We would therefore expect no offset between the spiral shock and star formation, and 
no obvious co-rotation radius. Spurs which occur in grand-design spirals when large 
clumps are sheared off leaving the spiral arms, are rare in the active, time-evolving 
spiral reported here. Instead, large branches are formed from spiral arms when the 
underlying spiral potential is dissolving due to the N-body dynamics. We find that the 
molecular cloud mass spectrum for the active potential is similar to that for clouds in 
grand design calculations, depending primarily on the ambient pressure rather than 
the nature of the potential. The largest molecular clouds occur when spiral arms 
collide, rather than by agglomeration within a spiral arm. 

Key words: galaxies: spiral - hydrodynamics - ISM: clouds - ISM: molecules - stars: 
formation - galaxies:structure 



1 INTRODUCTION 

■ A key objective to understanding star formation is to explain 
how molecular clouds and stars form in different types of 
, galaxies. The presence of spiral shocks in grand design galax- 
' ies explains the predominance of young stars in s piral arms 
; jRobert j Il969l: IwoodwardI Il976l : [Cepa fc BeckmanI Il990l : 



iKnapen et al.lll992l . 19961 ). and can also accou nt for the ob- 
served velocity dispersions in molecular clouds (IZhang et al.l 
I2OOII : iBonnell et Ed]l2006l : iKim. Kim fc Ostrikeijl200g i . Oth 
erwise, star forma tion is thought to occur through trigger- 
ing b y supernovae l|Mueller fc ArnettI 19761 : iGerola fc SeidenI 
1 19781 ) , spontaneously through gravitational, thermal or mag- 
ne tic instabilities, o r by turbulent compression (as reviewed 
in lElmegree"nl ll99gY 

Prom optical images, spiral galaxies can be classified as 
i) grand design, usually consisting of 2 symmetrical spiral 
arms, ii) multi-armed, with several asymmetric spiral arms 
or iii) flocculent, wit h multiple shorter arms. Using t he clas- 
sification scheme of lElmegreen fc Elmegre"enl (|l987h . grand 
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design galaxies are inclusive of multi-armed spiral galax- 
ies, but we retain this terminology to distinguish between 
the two types. Classical grand design structure is thought 
to originate from perturbations due to bars or companion 
galaxies which induce a density wave in the underlying stel- 
lar disk. Flocculent (and multi-armed) structure develops 
from local gravitational instabilities, which are sheared into 
short, transient spiral arms. Most fiocculent galaxies are 
also found to have a weak density wa ve in the older stel- 
lar population visible in the K band (jThornlev fc Mundvl 
1 19971 : lGrosb0l fc Patsisll 19981 ) . 

Several hydrodynamical simulations have in- 
vestigated the response of the ISM to spiral den- 
sity waves in grand design galaxies, by assuming 
a rigidly rotating spi ral potential, usually w it h an 
isothermal medium llPatsis. Grosbol fc Hiotelis 1997 : 



Chakrabarti. La ughlin fc Shul l2003l : iKim fc Ostrikeij bood : 
Dobbs fc Bonne ll 2006). The spiral pattern in these sim- 
ulations is assumed to be long-lasting, of at least several 
rotation periods. The formation of dense structures along 
the spiral arms are a s sociat e d with molecular clouds 
ijKim. Ostriker fc Stond l2003l : iDobbs. Bonnell fc Pringi3 
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I2OO6I : iDobbs fc Bonnelll l2007bh and the shearing of 
these features produc es spu rs perpendicular t o the 
arms dKim fc Ostrikerl l2002l : lOobbs fc Bonnelll l2006l : 



IShettv fc Qstrikerl |200{ ). For the purely hydrodyn amic 
simulations (IWada fc Kodal 12004 iDobbs et al] I2OO6I ). the 
development of molecular clouds and spurs is controlled 
by the dynamics of the shock. Simi lar features are formed 
through gravitational i nstabilities l|Kim fc OstrikeJ I2OO2I : 
IShettv fc Qstrikedl2006D . 

The generation of turbul ence has also been e xamined in 
simulations of galactic disks. IWada et al. I l|2002t ) show that 
gravitational instabilities combined with galactic rotation 
can drive turbulence, producing flocculent spiral structure. 
Super novae are expected to be a major progenitor of turbu- 
lence l|de Avillez fc Breitschwerdtl l2005l: I Joung fc Ma.c Low! 



I2OO6I ') in the ISM, whilst IPiontek fc Ostrikeij l|2005h sug- 



gest that MRI driven turbulence may be important in 
the outer parts of discs where supernovae are less fre- 
quent. An alternative possibility which has eme rged recently 



200ll:lBonnell et al.ll2006l: Kim et al.ll200e 


;lDobbs fc Bonnelll 


2007blV In particular. Bonnell et al.l (12006 


) show the passage 



of gas through a spiral arm and the formation of molecular 
clouds which exhibit a cr oc r"'^ velocity dispersion law. 

Realistic galactic potentials are not rigid, but time de- 
pendent. Even in grand design galaxies, the spiral pattern 
and pitch angle of t he spiral arms may change wit hin a 
few rotation periods l|Merrifield. Rand fc Meiddl2006f ). For 
galaxies where the spiral structure arises through more lo- 
calised gravitational instabilities, the number and shape of 
spiral arms are expected to be continually evolving. Numer- 
ous N-body simulations have investigated the spiral struc- 
ture arising from such instabilities in a stellar disc. If the 
mass of the halo is sufficient, the bar instability is pre- 
vented and the resulting structure then depends on the 
Toomre Q parameter. Generally gravitational instabilities 
produce flocculent and multi- armed spirals (,Bottcma 2003; 
iLi. Mac Low fc KlessenI l2005l) . unless an increase of Q in 
the centre of the disk is induced to suppress m > 2 modes 
iih omasson. Elmegreen. Donner fc Sundeliui Il990l ) . With 
Q ^ 1 — 2, multi-ar med spirals with modes up to m = 6 are 
formed l|Sellwood fc Car lbcrg Igsl). Lowering the disk mass 
or raising the veloci ty dispersion increases Q and p roduces 
a flocculent galaxy (|Elmegreen fc Thomassorilll993h . In all 
these cases, the gaseous component of the disc develops a 
much more pronounced spiral structure than present in the 
stellar component. Hydrodynamical simulations show that 
the gas becomes dense enough for star formation to occur 
in th e spiral arms, providing the Toomre criterion (iToomrd 
1 19641 ) is satisfi ed (Q < 1 for gas (|Tasker fc Brvanll2006l ) or 
gas and stars (|Li et al.ll2005l )). 

In this paper, we describe hydrodynamical simulations 
of multi-armed spiral galaxies, which use an active galac- 
tic potential. This p otential is t aken from an N- body sim- 
ulation described in ISellwood fc Carlbcrg] (|l984l ). The po- 
tential has also recentl y been used for 2D gas simulations 
l|Clarke fc Gittinsll2006D with the grid code CMHOG, which 
highlight the dependence of the star formation rate on the 
changing spiral structure. Whilst the gas dynamics in the 
grand design simulations are driven by a density wave, the 
potential from the N-body simulation arises through tran- 
sient gravitational instabilities in the disc stars. We describe 



the galactic structure, gas dynamics and the properties of 
molecular clouds which arise with the active potential, and 
compare with previous simulations wh ere a rigid potential 
was adopted l|Dobbs fc Bonnellll2007bl ). 



2 CALCULATIONS 

We use the 3D smoothed particle hydrodynamics (SPH) 
code based on the version by Benz (|Benz et al ]tl990'). The 
smoothing length is allowed to vary with space and time, 
with the constraint that the typical number of neighbours 
for each particle is kept near Nmigh ~ 50. Artificial viscosity 
is included with the standard parameters a = 1 and /3 = 2 
iMonaghan fc Lattan"ziolll985l : lMonaghanlll992l '). 

We describe the galactic potential and initial conditions 
for the simulations in Sections 2.1 and 2.2. The 3 simula- 
tions performed are also summarised in Table 1. In all the 
simulations, self-gravity magnetic fields, heating, cooling or 
feedback from star formation are not included. This paper 
focuses on the differences in the gas dynamics and structure 
resulting from the active potential. 



2.1 Galactic potential 

The gravitational potential for these calc ulations is taken 
from a previous N-body simulation ( Sellwood fc Carlberd 
1 1984 ). The N-body calculation includes both a disk and a 
rigid halo, the latter suppressing the formation of a bar in 
the disk. In the N-body simulation, the disk evolves to pro- 
duce an asymmetric spiral pattern with several spiral arms. 
The distribution of partic les in the calculation are sho wn at 
different times in Fig. 2 of lSellwood fc Carlberd l|l984 . The 
amplitude of the perturbation peaks after approximately 1 
rotation period, measured at the midpoint of the disk, after 
which the intensity of the spiral pattern quickly decays. The 
spiral perturbation to the potential decays in the N-body 
simulation as the velocity dispersion increases, and the disc 
becomes stable (i.e. Q>1) against gravitational instabilities. 

The gravitational potential and accelerations were de- 
rived fr om the N-body simulatio n and used in previous re- 
sults bv lClarke fc GittinsI (|2006l ). Both quantities were cal- 
culated over a 2D cylindrical polar grid consisting of 90 x 
128 points, after each time step in the N-body simulation 
(giving a total of 1500 grids) . In this paper we just use the 
accelerations, which we interpolate both spatially and with 
time. Thus at a given time in the SPH calculation, we use 
the accelerations from the polar grids at the closest preced- 
ing and succeeding times. 

Our simulations are 3D, so we also apply accelerations 
acting solely in the z direction assuming an axisymmetric 
logarithmic disk potential, i.e.: 



Ipdisk 



dz 



where 



(1) 



where Rc = 1 kpc, uo = 220 km s~^ and Zq— 0.7 is a measure 
of the disc scale height. The gravitational field is not fully 
consistent, since the z component is not related to the spiral 
potential. Thus the generatio n of turbulence by 'flapping' 
motions of the spiral shocks (|Kim et al.ll2006l ), is not able 
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Run 


No. of 


Mass of cold 


Mass of warm 


Total mass 


Average S 


Mass of H2 after 




particles 


/'inn fT^in flK/l \ 
(iUU Kj gas (Mq) 


(iU Kj gas (M0j 




(M0 pc ) 


zUU iVlyr v^Vl0 ) 


A 


27 million 


9 XIO* 






3.2 


1.5 xlO** 


B 


27 million 


6 xlO* 


6 xlO* 


1.2 xlO'' 


3.9 


3.9 xlO* 


C 


4 million 





1.2 xlO^ 


1.2 xlO^ 


3.9 





Table 1 . Table showing the different simulations described, and the proportions of warm and cold gas for each. 



to occur in these models, requiring either an enhanced z 
component of the potential in the spiral arms, or self gravity. 
A fully consistent potential may also increase the density in 
the spiral arms, again due to the increased axial component 
of the potential there. 



2.2 Initial distribution of gas 

The particles are initially distributed uniformly across a 
disk with a scale height of 100 pc and maximum radius 
of 10 kpc. We exclude particles from the inner 1 kpc ra- 
dius of the disk, but as the simulation progresses, particles 
may enter this region. The initial v elocities of the parti- 
cles as sume the rotation curve given in lSellwood fc Carlbergj 
||1984 ) (Equation 2), with a maximum rotational velocity of 
220 km s~^, and the particles are initially placed on circu- 
lar orbits. We choose a length scale such that the maximum 
extent of the potential is 1 4 kpc (corresponding to r=7 in 
ISellwood fc Carlberd l| 19841 11. although since the spiral per- 
turbation does not extend to the whole of the disk in the 
N body simulation, we adopt a maximum radius of 10 kpc. 
With this choice of length and velocity scales, the time for a 
simulation to complete the equivalent of the 1500 timesteps 
of the N-body calculation is 430 Myr. This corr e spond s to 
a time of t=4 on Fig. 2 of ISellwood fc Carlberd l|l984h (or 
equivalently 0.94 Gyr in lClarke fc Gittinj|2006l l! 

In addition to the rotational velocities, we also imposed 
a dispersion on the velocity distribution. The added disper- 
sion is selected from a Gaussian distribution with a standard 
deviation of 2.5% of the orbital speed. This was used for all 
the X, y, and z velocities. The velocity dispersion in the z 
direction maintains vertical equilibrium in the disk. Turbu- 
lence is not explicitly included in the initial conditions of 
our simulations. However, as mentioned in the introduction, 
previous results have shown that spiral shocks induce a sig- 
nificant velocity dispersion in the gas. Thus the turbulence 
is generated self consistently in the simulation as gas is sub- 
ject to spiral shocks, and the initial velocity dispersion is 
not very significant to the outcome of our results. 

We perform 2 high resolution simulations, with 27 mil- 
lion particles (Table 1), where the resolution is approxi- 
mately 40 Mq per particle. One calculation contains pre- 
dominantly cold gas, with a much smaller fraction of warm 
gas, whilst the other uses equal phases of warm (lO"* K) and 
cold gas. The SPH particles are initially randomly assigned 
as warm or cold gas. We also perform a lower resolution cal- 
culation with just warm gas, and in all calculations the gas 
is isothermal. 



2.3 Two phase medium 

In this present study, we neglect any heating or cooling be- 
tween the 2 phases (so cold particles remain cold through- 
out the simulation). Although a full thermal treatment is 
preferable, a two phase medium is a first approximation 
to modeling cold and warm gas simultaneously, and al- 
lows a direct comparison with our previous calculations 
l|Dobbs fc Bonnelll l2007bl l. Observations indicate both of 
these components are present in the ISM (|Heiles fc Trolandl 
l2003l l. and thus it is important to include them in spi- 
ral shocks. We assume that cold gas is pre-existing in the 
ISM (rather than the o utcome of cooling from the wa rmer 
phase), as suggested bv lPringle. Allen fc Lubowl l|200ll l. Al- 
though if the cooling timescale of w arm gas in the ISM 



ver y short (i.e. a few Myr, e.g. Audit fc Hennebelld 



2005; Heitsch ct al. 2006|; IVazauez-Semadeni et al.l l2007l : 
Glover fc Mac Low 2007i l cold clumps would quickly develop 



as warm gas enters the spiral arms, producing conditions in 
the spiral shock similar to these calculations. 

Although randomly positioned initially, the warm and 
cold phases become separated during the simulation. This 
separation is expected when gas passes throu gh a spiral 
shock. In p revious simulations (Dobbs fc Bonnelll l2007bl : 
iDobbd I2OO8I I. cold gas is compressed into well resolved 
clumps of order 10® M©. However phase separation is also 
apparent before the gas has time to pass through the shock. 
This is particularly noticeable at large radii in the calcula- 
tions presented here, where the dynamical time is greatest. 
This is a numerical rather than physical artifact, and oc- 
curs due to a discontinui ty in pressure between the warm 
and cold gas l|Price|[2007l l. The separation of cold gas into 
clumps is accentuated by the shear velocities of the gas. 
When cold particles are adjacent to each other, a pressure 
imbalance between the two phases causes the cold gas par- 
ticles to move closer together rather than move freely past 
each other. This pressure imbalance efi'ectively suppresses 
mixing between the two phases of gas. 

Apart from in the inner regions of the disk, gas flowing 
into the spiral arms is situated in cold clumps surrounded 
by a warm medium. This distribution is nevertheless much 
more similar to the physical conditions in the ISM, where 
the cold gas is situated within higher density clumps, rather 
than randomly interspersed with the warm gas. However the 
distribution and spacing between these clumps is dependent 
on the resolution, with clumps typically containing 50-100 
particles. Ideally we would wish to remove this numerical in- 
stability, ensuring clumps produced in the calculations are 
solely related to the spiral shocks. This should be feasible 
for future calcul ations, by w ay of an artificial thermal con- 
ductivity term |Prici HqoJ), which smoothes the pressure 
between adjacent gas particles of different temperatures. 
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2.4 Calculation of molecular gas density and 
cooling of molecular gas 

We estimate the fraction of molecular gas usi ng tlie equa- 
tion fo r the rate of change of H2 density from iBergin et alj 

(Eooi), 

^^^^ = Rgr{T)npn{H) ~ [Ccr + U.s{N{H2), Av)]n(H2). 

(2) 

The total number density is Up, where Up = n{H) + 2n(H2), 
N is the column density (of atomic or molecular hydrogen) 
and T the temperature. The photodissociation rate is 

Cdiss=CdiU0)f{NiH2),Av), (3) 

from iDraine fc Bertoldil (|l996l ). The constant Cdiss{0) = 
4.17 X lO"^'^ s~^ is the photodissociation rate for unshielded 
H2. The H2 is assumed to be subject to a radiation field 
equivalent to the U.V. radiation emitted by BO stars av- 
eraged over the plane of the Galactic disk. The function 
f{N{H2),Av) takes into account self shielding of the H2. 
This function includes the dust extinction (Ay), and the 
column density of the H2, the latter calculated using the 
local density of H2 and the scale height of the disk. The 
second destructive term is the cosmic ray ionization rate, 
(^cr = 6 X 10"^* s~^. The rate of formation on grains is 
Rgr, assuming an efficiency of S — 0.3. We do not explic- 
itly change S for the warm gas in our calculations, but the 
density of the warm gas is too low to prevent immed iate 
photodissociation of any H2 formed (jPobbs et al.ll200(l ). 

Although we do not include heating or cooling between 
the 2 phases, we do allow cooling of the 100 K gas which 
becomes molec ular. We use a simple polytropic equation of 
state P oc p° '' (|Larsonll2 005') suggested for molecular clouds. 
We take only the molecular gas density for this equation 
of state, setting a maximum temperature of 100 K for gas 
which is fully atomic and a minimum temperature of 30 
K for gas which is fully molecular. It is not thought that 
the addition of molecular cooling has a significant effect on 
these simulations, but may increase the density slightly in 
the shocks and prolong the time gas remains molecular. As 
the molecular gas dissociates, the gas heats back up (to a 
maximum temperature of 100 K) using the same equation 
of state. 



3 RESULTS 

The primary driver for the evolution of the gas is the devel- 
opment of spiral arms in the N-body potential. The main 
difference between the simulations presented here and those 
for a grand design galaxy is that the potential changes on 
timescales comparable to the gas crossing time of a spiral 
arm. Rather than gas passing through a spiral shock, the 
evolution and structure of the gas is determined by local 
changes in the potential minimum and collisions between 
spiral arms. 

3.1 Structure of the disc 

The development of the spiral structure with time for the 3 
different calculations is shown in Fig. 1. The perturbation 
to the potential, and hence the extent of the spiral arms, is 



at first restricted to the inner disc. The spiral arms spread 
to larger radii as the magnitude and extent of the spiral 
perturbation to the potential increase. The number of spiral 
arms is not constant and the dominant mod e varies between 
m = 2 and m = 8 (| Clarke fc GittinsI 1200611. The structure 
also changes as spiral arms collide with each other, giving 
the appearance of bifurcations of the arms into different 
branches. Enhanced star formation is pre dicted where high 
densi ty regions develop as the arms collide (|Clarke fc Gitting 
I2OO6I ). Approximately half way through the N-body simula- 
tion, the spiral perturbation in the disk begins to decrease. 
Consequently the spiral structure begins to wind up after 
about 200 Myr in our simulations. 

The large scale structure and spiral features (Fig. 1) 
are similar for the different thermal distributions. With cold 
gas present there are dense clumps in the shocks and higher 
density regions, whereas when there is only warm gas (lower 
panels), the spiral arms are smooth. In some regions, the 
gas appears to form 'double shocks', much more noticeable 
for the cold gas, altho ugh still apparent for th e simulation 
with warm gas and in IClarke fc Gitting l|2006l ) . These are 
due to new spiral arms forming, which attain similar rota- 
tional velocities to older spiral arms. The structure towards 
the centre of the predominantly cold simulation becomes 
chaotic, due to repeated collisions of stronger shocks in this 
region. 

A detailed section of the simulation with 50% warm gas 
(Run B) is shown in Fig. 2. There is considerable structure 
in the cold gas in both the spiral arms and inter-arm re- 
gions. However, the agglomeration of gas into larger clumps 
in the shock is not as apparent as when ther e is a standin g 
density wave from a grand design potential l|Dobbdl2008l ). 
Instead, collisions between spiral arms are more important 
(e.g. Fig 2), and this is where the largest structures form in 
the gas (and thus where agglomeration of clumps is occur- 
ring). 

3.2 Gas dynamics of dynamically evolving galaxy 

The response of the gas to the active potential appears to 
be quite different to the behaviour of gas subject to a grand 
design potential. We describe here the main differences be- 
tween the two cases. 

For the active potential considered here, the gas essen- 
tially traces the potential minimum. This is in contrast to 
grand design galaxies, where the gas shocks as it passes the 
minimum but then leaves the spiral shock to continue to the 
next spiral arm. Here gas enters the minima of the potential 
as they form, unlike the grand design case where the min- 
ima already exist. The gas then tends to accumulate in the 
potential minima rather than climbing out of the potential, 
retaining spiral structure even as the potential minimum 
dissolves. 

In Fig. 3, we show the typical evolution of gas as it 
enters a potential minimum. The figure highlights a selection 
of particles at 4 different time frames from the simulation 
with warm gas. The particles were selected as those in the 
section of spiral arm shown in the third frame, with their 
locations at earlier times shown in the top two panels. As 
a minimum develops in the underlying potential, a shock 
occurs as gas falls into the minimum and produces a spiral 
arm. The top two panels of Fig. 3 indicate that particles 
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Figure 1. The time evolution of the spiral pattern is shown in these column density plots, which are inclusive of both cold and warm 
phases of gas. The top row shows Run A with predominantly cold gas, the middle with equal proportions of cold and warm gas (Run B) 
and for the bottom row, all the gas is 10'* K (Run C). Each plot shows a region of -11 kpc< x <llkpc and -11 kpc< y <llkpc, assuming 
a Cartesian axis centred at the origin, and gas is flowing anticlockwise round the disk. The main spiral arm features are analogous for 
the different thermal distributions, although there is much more structure on smaller scales when there is cold gas. 
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Figure 2. These images zoom in on a small region from tlie simulation witli 50% cold gas, Run B. The plots correspond to a time of 200 
Myr, and show areas of size i) 22 kpc x 22 kpc, ii) 8 kpc x 8 kpc, iii) 4 kpc x 4 kpc and iv) 1.5 kpc x 1.5 kpc. The gas appears increasingly 
structured at smaller scales, becoming sheared into filamentary structures in the inner regions. The figure focuses on a bridge of gas 
between two spiral arms which corresponds with the remnants of a previous arm. Dense material has agglomerated either side of this 
bridge, where spiral arms have collided. The same scale is used as shown in Fig. 1. 



enter the minimum from both sides of the developing spiral 
arm. The spiral arm in the bottom left panel is much sharper 
on the upst ream side where the gas generally flows into the 
spiral arm (iBonnell et al.l[2006l ). The fourth panel in Fig. 3 
(bottom right) shows the particles at a much later time, by 
which time they have formed a completely different spiral 
arm. 

Fig. 4 shows the densest gas from 2 of the simulations 
overplotted on the underlying potential. In the top panel, 
the gas is taken from the simulation with predominantly 
cold gas (Run A) and for the other, the gas is warm (Run 
C). Only the non-axisymmetric component of the poten- 
tial is shown, calculated at each point by subtracting the 
mean value of the potential over the corresponding cylin- 
drical radius. In both cases the high density gas strongly 
coincides with the potential minimum. Since gas effectively 



falls into the spiral potential, this is not surprising. The co- 
incidence between the spiral arms and the potential min- 
ima is in contrast with grand design models where the spi- 
ral shoc k s ten d to be offset from the potential minimum. 
iRobertj (| 19691 ) found that spiral shocks occur on the lead- 
ing side of the potential minimum inside co-rotation for 
quasi-stationary solutions with m = 2. However the location 
varies depending on sound speed, pattern speed, number 
of spiral arms and magnetic field str ength, as shown in re- 
cent simulations _^yz, Kranz & Rix 20031: ICittins fc Clarke! 
l2004l : iDobbs fc Pric6.200& ,Wada.2007i ). For the active po- 
tential applied here, a quasi-stationary state of course does 
not evolve, since the spiral potential changes with time, and 
we would not expect an offset between the spiral shock and 
star formation. 

There are nevertheless some sections of spiral arms that 
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Figure 3. The particles are shown for a subsection of the disc at times of 70 (top left), 76 (top right), 82 (bottom left) and 195 (bottom 
right) Myr for the simulation with warm gas (Run C). Particles indicated in black were selected from the spiral arm displayed in the 
bottom left panel, and are shown at earlier and later times in the other panels. The first three panels show the accumulation of gas into 
a spiral arm. The spiral arm structure of the gas persists as the gas encounters multiple spiral arm passages, as indicated by the fourth 
panel. At the latest time, only 1 in 10 of the total number of particles are plotted. 



appear to cross between two minima.The stellar spiral arms 
associated with the N-body simulation are transient, so the 
minimum will eventually disappear. However the gaseous 
spiral arms still retain their spiral structure, even though 
the minimum dissolves. Eventually the gas reaches another 
minimum and joins a new spiral arm. The final panel in 
Fig. 3 shows a different section of spiral arm at a much 
later time frame. The selected particles have already passed 
3 spiral arms at this point. Gas particles highlighted with 
y > —2 kpc are largely coincident with a potential mini- 
mum, and also corresponds to one of the dense arms of gas 
in Fig. 4. On the other hand, gas particles with x > —2 kpc 
lie between two minima or are forming the new spiral arm 



seen in the lower part of the figure. In Fig. 4, spiral arms 
lying between the potential minima are most evident when 
cold gas is present and the densities of the gas are higher. 
A further consequence of gas retaining the spiral structure 
imposed by the potential is that at later times in the simu- 
lations (e.g. 3rd time frame. Fig. 1) , most of the gas is in the 
spiral arms. Regions between the arms are relatively empty, 
and since little gas is entering the potential, the spiral arms 
correspond merely to dense regions of gas rather than actual 
spiral shocks. 

In addition to the density, the velocity of the gas which 
is in the spiral arms largely follows that of the stellar com- 
ponent of the disc, with different pattern speeds emerging 
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Figure 4. The most dense gas from simulations with predomi- 
nantly cold (Run A, top) and warm (Run C, lower) is overplotted 
on the potential from the N-body simulation. Only the spiral 
component of the potential is indicated, and black/dark red indi- 
cate regions where the potential is lowest, yellow/ white indicates 
the potential is highest. The gas is shown in blue. Generally the 
densest gas corresponds with the minima of the potential. The 
regions shown are 20 kpc x 20 kpc and the potential rotates in 
the anti-clockwise direction. The corresponding time is 200 Myr. 



for different sections of spiral arm. The rotational velocities 
of the gas thus indicate that there is no fixed pattern speed 
or co-rotation radius for the spiral perturbation. 



3.3 Velocity dispersion 

Fig. 5 shows the velocity dispersion for the 2 phase simu- 
lation with equal components of warm and cold gas (Run 
B) as a function of azimuth. The average velocity dispersion 
is calculated over a ring of width 200 pc and azimuthally 
divided into 100 segments. The ring is centred at 2 different 
radii, 4 kpc and 8 kpc. Both show peaks where the gas passes 
through spiral shocks, but whereas at radii of 4 kpc, the ve- 
locity dispersion reaches > 10 km s~^, at 8 kpc the velocity 
dispersion barely exceeds 3 km s~^. All the gas (warm and 
cold) is used to determine the velocity dispersion in Fig. 5. 



0) 

> 




Figure 5. The velocity dispersion is shown against azimuth after 
200 Myr for Run B (50% cold gas). The angle d is measured 
anticlockwise, in the direction of the flow, at radii of 4 and 8 kpc. 



The variation in the velocity dispersion is dominated by the 
cold gas. 

The increase in the velocity di spersion occurs as clumpy 
gas passes through a spiral shock (|Dobbs fc Bonnellll2007ah 
and the magnitude of the increase is dependent on the Mach 
number of the shock. At smaller radii, the potential is much 
stronger, so the Mach number and therefore velocity dis- 
persion are higher. Likewise, for a given radius, the velocity 
dispersion in the spiral arms increases with time whilst the 
strength of the perturbation to the potential increases. How- 
ever the velocity dispersion starts decreasing by 250 Myr, 
since little gas is entering the spiral arms and shocks are no 
longer occurring. 



3.4 Spurs and branches 

In simulations of grand design galaxies, regular spurs are 
found to extend from the arms into the inter-arm regions. 
These spurs are due to the shearing of clumps of gas leaving 
the spiral arms. They ca n occur if there is cold g as (< 1000 
K) without self gravity (jPobbs fc Bonnelill2006l), or evolve 
from gravitational instabilities in warm gas ( Kim fc Ostriked 
I2OO6I ). These features are not apparent in the simulations 
with the active potential. As mentioned in Section 3.2, gas 
tends to stay in the spiral arms and retain spiral structure, 
even as the minimum dissolves. Thus without gas leaving the 
spiral arms, gas cann ot form spurs. From a survey of feather- 
ing in spiral galaxies (|La Vigne et al.ll2006l ). 20% of galaxies 
are found to have periodic feathers. Generally galaxies with 
complex or fiocculent structures are less likely to contain 
clearly defined feathers, although those with clearly defined 
dust lanes are more likely to contain such features. 

Instead, when a minimum in the potential disappears, 
the gas retains the shape of the whole arm. Thus much larger 
branches occur in these simulations, located between the 
main spiral arms. A clear example is shown in Fig. 2b. These 
features are not present in simulations with ti me indepen- 
dent potentials, except where r esonances occur (jPatsis et al.l 
I1997I : IChakrabarti et ahl |2003| ) . As the potential minimum 
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Figure 6. These images display tlie column density of molecular hydrogen (red) overplotted on the total column density of HI, from 
Runs A and B. The left panels show Run A, with 90% cold gas, and the right Run B, with 50% cold gas. The top panels show regions 
of 16 kpc X 16 kpc, and the bottom panels show an area of 4 kpc x 4 kpc (where -3 kpc< x <1 kpc and 2 kpc < y <6 kpc assuming a 
Cartesian grid) the same as in Fig. 2. All the panels correspond to a time of 200 Myr. In both simulations the molecular gas is situated in 
highly structured clumps. With a larger component of warm gas, there is more inter-arm molecular gas, which is sheared into filamentary 
clumps. 



decays, the velocity of the spiral arm adjusts to that of cir- 
cular velocities, and the shape of the spiral arm changes by 
differential rotation. Large branches can be seen extending 
from the main spiral arms in many galaxies such as MlOl, 
M74 and M33. The pitch angle and length of these features 
tend to vary much more than would be expected from spurs 
which originate from instabilities in the spiral arms. 



3.5 Molecular gas in the disk 



We show the distribution of molecular gas, overplotted 
on the total gas density in Fig. 6, for Runs A and B. 
Similar to previous resu lts with a grand design potential 
dPobbs fc Bonnellll2007bl ). there is much more molecular gas 
when there are equal amounts of cold and warm gas, com- 
pared to when the gas is nearly all cold. This is because, 
as mentioned in the previous paper, more of the inter-arm 
gas is molecular, when it is confined to higher densities by 
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the warm gas. There is also more molecular gas in these 
simulations compared to the grand design case. This is due 
to a combination of factors: a stronger potential, the fact 
that gas tends to remain in the spiral arms and collisions 
between spiral arms leading to more gas at higher densities. 
After 200 Myr (Table 1), 40% and 23% of the total disk mass 
is molecular for Runs A and B. However the system is not 
yet in equlibrium, and both fractions are still increasing. 

3.6 Molecular clouds 

We also applied a clump-finding algorithm to the distribu- 
tion of molecular gas in Runs A and B, predominantly to 
calculate the mass spectrum for comparison with previou s 
results for grand design galaxies l|Dobbs fc Bonnellll2007bh . 
For Run B, from Fig. 6, it is apparent that clumps in the 
outer part of the disk, where there is no spiral perturbation, 
have formed which are numerical in origin. At later times 
in the simulations, these clumps have retained high densi- 
ties for a sufficient time that they contain molecular gas. We 
therefore only consider a region of 10 kpc by 10 kpc centred 
at the midpoint of the disk, which does not include these 
clumps. 

As for iDobbs fc Bonnelll (|2007bl ). we use a clump- 
finding algorithm to identify molecular clouds. We take 
the clump-fin ding method CFl previously described in 
iDobbs fc Bon ncU (2007b!}, which selects grid cells above a 
certain column density and classes groups of adjacent cells 
as a single clump. We select cells which have a molecular 
gas surface de nsity above 4 Mp pc as used for our pre- 
vious results l|Dobbs fc Bonnelll l2007bl ) , and only include 
clouds with > 30 particles. The resolution of the simula- 
tions shown here is ove r twice that of previous analysis 
l|Dobbs fc Bonnellll2007bh with the mass of each particle ap- 
proximately 40 Mq. The minimum total mass of a cloud is 
then 10'' M0, although the mass of molecular gas may be 
less. 

We find that the most massive molecular clouds (> 
10^ Mq) are located where spiral arms are colliding or 
merging. This agrees with the previous predictions by 
IClarke &: Gittina (|200f?) that enhanced star formation would 
occur here. We also find that there are very few inter-arm 
clouds between 2 and 4 kpc where most of the gas is in the 
spiral arms, whereas at lower radii there are more interarm 
clouds as a result of multiple collisions between spiral arms. 

Fig. 7 shows the mass spectra for Runs A and 
B. Assuming a mass spectrum of dN/dM — M~°', 



a — 1.75 ± 0.12 for the simulation with predominantly cold 
gas, and a = 2.05 ± 0.1 with a larger (50%) component 
of hot gas. The mass distribution aris es through cloud 



coaeulation (e.s. Field & Saslaw 


19651: Taff & SavedoflS 


19721: iHandburv et al. 


Il977l: 


ScoviUe & HershI 119791: 


Kwan & Valded 19871: 


Das & Joe 


19961: Clark & Bonnelll 


20061: iHarfst et all l20od 


), a process enhanced by collisions 



between spiral a rms. For previous simula tions of grand 
design galaxies (|Dobbs fc Bonnellj l2007lj ). we obtained 
similar values and also found that the mass spectrum for 
molecular clouds was steeper with a larger component of 
warm gas. This suggests that the nature of the ISM is more 
important in determining the mass spectrum than the na- 
ture of the potential, and how clump agglomeration occurs. 
Previous observations of the Milky Way found a = 1.5 — 1.8 




Run A - 
Run B l 



0-5 r nHj- h 

0.0 F, 1, 1, 1, 1,1,1, 1, 1, III], 1, 1, 1, 1, 1,1,1, III, 1, 1, 1, 1,1,0, 1, 1, Ij.D.Oi , ,U - 

2 3 4 5 6 7 

Log Mass (M^) 

Figure 7. The clump mass spectra are shown for Runs A and B. 
Both are determined using the molecular mass of the clumps, after 
200 Myr. The slopes indicated correspond to a = 1.75 (Run A) 
and a = 2.05 (Run B) where a is the index of the mass spectrum. 
The mass spectrum is steeper when more warm gas is present. 

(jSolomon et al.l Il987f ) , although some more recent survey s 
of external galaxies (iBlitz et all l2007l : iLada et all |2007| ) 
indicate that a may be higher. 



4 CONCLUSION 

We have extended previous analysis of the structure and 
dynamics of the ISM in spiral galaxies by applying a time- 
dependent spiral potential derived from an N-body calcu- 
lation. The gas forms long spiral arms, exhibiting a multi- 
armed structure where the underlying stellar component is 
determined by gravitational instabilities rather than a spiral 
density wave. The main difference compared to the grand 
design case is that the development, collision and dissipa- 
tion of spiral arms in the underlying potential is the main 
driver for generating high density structure in the ISM. The 
timescale for the evolution of structure is thus linked to the 
lifetime of the local minimum. 

Gas accumulates in the potential minima as they de- 
velop, thus the densest gas is coincident with the stellar 
spiral arms. Consequently no ofi^set would be expected be- 
tween the dust lanes and star formation in such a galaxy. 
Spurs associated with the shearing of GMCs as they leave 
the spiral arm, are largely absent with the active potential. 
This is because gas only leaves the spiral arms once the po- 
tential minimum dissipates. The gas then retains the large 
scale spiral arm structure, becoming slowly wound up until 
it fully dissipates or encounters another spiral arm. Instead 
of spurs, much larger bridges, which are remnants of former 
spiral arms, extend between the main spiral arms. 

We further identified molecular clouds from these sim- 
ulations, finding that the largest GMCs occur where spiral 
arms are colliding. We found the surprising result that the 
mass spectra are similar to previous results where the po- 
tential corresponded to a density wave. This suggests that 
the external pressure applied to the cold gas determines the 
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sizes and masses of the clouds, with a steeper index when 
the cold clumps are pressure confined. 

A natural extension of this work would be to model both 
the stellar and gaseous components of the disc, to include 
the gravitational instabiliti es of the stars and the response 
of the gas simultaneously l|Li et al. I [2005.). Furthermore in- 
corporating stellar feedback is required to treat the stellar 
and gaseous components self consistently. We have also ne- 
glected heating and cooling, magnetic fields and self gravity, 
which we are in the process of analysing with respect to 
grand design galaxies. We leave the consideration of these 
processes for more general stellar potentials to future work. 



ACKNOWLEDGMENTS 

We are grateful to Jerry Sellwood and Dave Gittins for pro- 
viding data from the N-body simulation described in Sec- 
tion 2.1. We also thank the referee. Ant Whitworth, for help- 
ful comments which significantly improved the clarity of the 
paper. Computations included in this paper were performed 
using the UK Astrophysical Fluids Facility (UKAFF). 

This work, conducted as part of the award 'The for- 
mation of stars and planets: Radiation hydrodynamical and 
magnetohydrodynamical simulations' made under the Eu- 
ropean Heads of Research Councils and European Science 
Foundation EURYI (European Young Investigator) Awards 
scheme, was supported by funds from the Participating Or- 
ganisations of EURYI and the EC Sixth Framework Pro- 
gramme. 



REFERENCES 

Audit E., Hennebelle P., 2005, A&A, 433, 1 

Benz W., Cameron A. G. W., Press W. H., Bowers R. L., 

1990, ApJ, 348, 647 
Bergin E. A., Hartmann L. W., Raymond J. C, Ballesteros- 

Paredes J., 2004, ApJ, 612, 921 
Blitz L., Fukui Y., Kawamura A., Leroy A., Mizuno N., 

Rosolowsky E., 2007, in Reipurth B., Jewitt D., Keil K., 

eds, Protostars and Planets V, Univ. Arizona Press, Tus- 

con, p. 81 

Bonnell I. A., Dobbs C. L., Robitaille T. P., Pringle J. E., 

2006, MNRAS, 365, 37 
Bottema R., 2003, MNRAS, 344, 358 
Cepa J., Beckman J. E., 1990, ApJ, 349, 497 
Chakrabarti S., Laughlin C, Shu F. H., 2003, ApJ, 596, 

220 

Clark P. C, Bonnell I. A., 2006, MNRAS, 368, 1787 
Clarke C, Gittins D., 2006, MNRAS, 371, 530 
Das M., Jog C. J., 1996, ApJ, 462, 309 
de Avillez M. A., Breitschwerdt D., 2005, A&A, 436, 585 
Dobbs C. L., 2008, MNRAS, submitted 
Dobbs C. L., BonneU I. A., 2006, MNRAS, 367, 873 
Dobbs C. L., BonneU I. A., 2007a, MNRAS, 374, 1115 
Dobbs C. L., BonneU I. A., 2007b, MNRAS, 376, 1747 
Dobbs C. L., Bonnell I. A., Pringle J. E., 2006, MNRAS, 
371, 1663 

Dobbs C. L., Price D. J., 2008, MNRAS, 383, 497 
Draine B. T., Bertoldi F., 1996, ApJ, 468, 269 



Elmegreen B. G., 1996, in Blitz L., Teuben P. J., eds, lAU 
Symp. 169: Unsolved Problems of the Milky Way, Kluwer, 
Dordrecht, p. 551 
Elmegreen B. G., Thomasson M., 1993, A&A, 272, 37 
Elmegreen D. M., Elmegreen B. G., 1987, ApJ, 314, 3 
Field G. B., Saslaw W. C, 1965, ApJ, 142, 568 
Gerola H., Seiden P. E., 1978, ApJ, 223, 129 
Gittins D. M., Clarke C. J., 2004, MNRAS, 349, 909 
Glover S. C. O., Mac Low M.-M., 2007, ApJ, 659, 1317 
Grosb0l P. J., Patsis P. A., 1998, A&A, 336, 840 
Handbury M. J., Simons S., Williams I. P., 1977, A&A, 
61, 443 

Harfst S., Theis C, Hensler G., 2006, A&A, 449, 509 

Heiles C, Troland T. H., 2003, ApJ, 586, 1067 

Heitsch F., Slyz A. D., Devriendt J. E. G., Hartmann L. W., 

Burkert A., 2006, ApJ, 648, 1052 
Joung M. K. R., Mac Low M.-M., 2006, ApJ, 653, 1266 
Kim C.-G., Kim W.-T., Ostriker E. C, 2006, ApJL, 649, 

L13 

Kim W., Ostriker E. C, 2002, ApJ, 570, 132 
Kim W.-T., Ostriker E. C, 2006, ApJ, 646, 213 
Kim W.-T., Ostriker E. C, Stone J. M., 2003, ApJ, 599, 
1157 

Knapen J. H., Beckman J. E., Cepa J., Nakai N., 1996, 
A&A, 308, 27 

Knapen J. H., Beckman J. E., Cepa J., van der Hulst T., 

Rand R. J., 1992, ApJL, 385, L37 
Kwan J., Valdes F., 1987, ApJ, 315, 92 
La Vigne M. A., Vogel S. N., Ostriker E. C, 2006, ApJ, 

650, 818 

Lada C. J., Alves J. F., Lombardi M., 2007, in Reipurth B., 

Jewitt D., Keil K., eds, Protostars and Planets V, Univ. 

Arizona Press, Tuscon, p. 3 
Larson R. B., 2005, in CorbeUi E., Palla F., Zinnecker H., 

eds, ASSL Vol. 327: The Initial Mass Function 50 Years 

Later, Springer, Dordrecht, p. 329 
Li Y., Mac Low M.-M., Klessen R. S., 2005, ApJ, 626, 823 
Merrifield M. R., Rand R. J., Meidt S. E., 2006, MNRAS, 

366, L17 

Monaghan J. J., 1992, ARA&A, 30, 543 
Monaghan J. J., Lattanzio J. C, 1985, A&A, 149, 135 
Mueller M. W., Arnett W. D., 1976, ApJ, 210, 670 
Patsis P. A., Grosbol P., Hiotelis N., 1997, A&A, 323, 762 
Piontek R. A., Ostriker E. C, 2005, ApJ, 629, 849 
Price D. J., 2007, ArXiv:0709.2772 

Pringle J. E., Allen R. J., Lubow S. H., 2001, MNRAS, 
327, 663 

Roberts W. W., 1969, ApJ, 158, 123 

ScoviUe N. Z., Hersh K., 1979, ApJ, 229, 578 

Sellwood J. A., Carlberg R. G., 1984, ApJ, 282, 61 

Shetty R., Ostriker E. C, 2006, ApJ, 647, 997 

Slyz A. D., Kranz T., Rix H.-W., 2003, MNRAS, 346, 1162 

Solomon P. M., Rivolo A. R., Barrett J., Yahil A., 1987, 

ApJ, 319, 730 
Taff L. G., Savedoff M. P., 1972, MNRAS, 160, 89 
Tasker E. J., Bryan G. L., 2006, ApJ, 641, 878 
Thomasson M., Elmegreen B. G., Donner K. J., Sundelius 

B., 1990, ApJL, 356, L9 
Thornley M. D., Mundy L. G., 1997, ApJ, 490, 682 
Toomre A., 1964, ApJ, 139, 1217 

Vazquez-Semadeni E., Gomez G. C, Jappsen A. K., 
Ballesteros-Paredes J., Gonzalez R. F., Klessen R. S., 



© 0000 RAS, MNRAS OOO.ITM 



12 C. L. Dobbs, I. A. Bonnell 



2007, ApJ, 657, 870 
Wada K., 2007, ArXiv:0710.0407 
Wada K., Koda J., 2004, MNRAS, 349, 270 
Wada K., Meurer G., Norman C. A., 2002, ApJ, 577, 197 
Woodward P. R., 1976, ApJ, 207, 484 
Zhang X., Lee Y., Bolatto A., Stark A. A., 2001, ApJ, 553, 

274 

This paper has been typeset from a TiTjX/ ffl^rjX file prepared 
by the author. 



© 0000 RAS, MNRAS OOP. [HO 



